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ABSTRACT 

I derive a set of canonical elements that are useful for collision orbits (peri- 
helion distance approaching zero at fixed semimajor axis). The coordinates are 
the mean anomaly and the two spherical polar angles at aphelion. 

The purpose of this short note is to derive a set of canonical orbital elements that are well- 
behaved for highly eccentric Kepler orbits. More precisely, I wish to find a set of elements 
that remain well-defined for collision orbits (eccentricity e — > 1 at fixed semimajor axis a). 
Such elements may be particularly useful for describing the orbital evolution of long-period 
comets. 

The Delaunay elements are defined as 

L={M*a) 1/2 , G = L(l-e 2 ) 1/2 , H = Gcosi, £, g, h; (1) 

here is the mass of the star in units where the gravitational constant is unity; a, e, and i 
are the semimajor axis, eccentricity, and inclination; and £, g and h are the mean anomaly, 
argument of perihelion, and longitude of the ascending node. The vector G is the specific 
angular momentum, G = |G| is its magnitude, and H = G • z is its z-component. In these 
variables the Kepler Hamiltonian is 

M 2 

Hk(L) = (2) 



Any set of six scalar orbital elements is indeterminate for some orbits (e.g. [Stiefel & 
|Scheifele 1971| ). For example, Delaunay elements are unsuitable for circular orbits [I and 
g indeterminate), equatorial orbits (g and h indeterminate), and collision orbits (i, g and 
h indeterminate). Poincare elements are suitable for circular orbits; eccentric and oblique 
elements are suitable for circular and for equatorial orbits; but all of these element sets 
are indeterminate for collision orbits (Murray & Dermott 1999). Herrick (1971) gives an 
extensive discussion of universal variables, which are well-defined for all Keplerian orbits 
including collision orbits, but Herrick's variables are not canonical. 
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The orientation of the line of apsides of an orbit can be specified by the spherical polar 
coordinates (6 a , <fi a ) of the unit vector a pointing from the origin to aphelion, which are given 
in terms of the elements i, g and h by 

cos# a = — sing sin sin 6 a sin(0 a — h) = — sing cos sin^a cos(0 a — h) = — cosg. (3) 

The angles (8 a , (f) a ) remain well-defined for collision orbits and thus are attractive candidates 
for canonical coordinates. I now demonstrate that the following is a set of canonical elements: 

l, e, h, e, e a , <p a , (4) 

where G is the component of the angular momentum along a line perpendicular to the z-axis 
and the line of apsides. More specifically, if a is the unit vector pointing towards aphelion, 
and we define a unit vector t = z x &/ sm9 a = (— sin0 a , cos0 a , 0), then = G • t, thus 
specifying the sign of O. The elements 9 a and (f) a are indeterminate for circular orbits, and 
is indeterminate if the line of apsides coincides with the z-axis. 

Consider a canonical transformation from the Delaunay elements (]T|) to new momenta 
and coordinates (I, w), defined by the mixed- variable generating function 

S(L,G,H,w) = wxL + (w 3 - \it)H - \kG 

_i Gcosu> 2 u -i Hcotw 2 

±GC ° S (G 2_ff 2) l/2 ^ HCOS (G 2_ff 2) l/ 2 > ( 5 ) 

where < W2 < vr and < W3 < 2n. The arguments of the inverse trigonometric functions 
have absolute value less than unity if 

sinw 2 > I cosi|, (6) 
that is, u>2 lies between ^ir — i and \ti + i. Then 

= , 2 ^. ± ( G2 _J^) V2 , I,™* (7, 

owi ow 2 \ sin W2 J 0W3 



and 



OS OS , Gcosw 2 

— Wi, g = — - = —^n ± cos 



dL x ' y 8G 2 (G 2 -H 2 y/ 2 ' 

dS 1 _j Hcotw 2 . . 

h = dH = W3 ~ 27lTCOS (G 2 -H 2 y/ 2 ' () 

Thus —n < g + \ix < 7r, and the upper sign applies when g + |7r > and the lower when 
g+ \ix < 0. 
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The equations for g and h imply that 



sin(u> 3 — h) = coticotw 2 



sin g cos i 



cos W2 = — sin g sin 2, 



sin W2 



cos(tf3 — h) 



=F(1 — cot 2 % cot 2 W2) 



1/2 



cosg 



(9) 



sin w 2 



Comparison with equations (§) then shows that w 2 = a and w 3 = <p a . It is now straightfor- 
ward to show that J 2 = 9. 
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